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I. INTRODUCTION 
, A. Generalized parton distributions 

o . 
o ■ 

■ The theoretical description of strong interactions in hard processes is based on the concept of QCD factorization. In 
many hard processes the quark-gluon interaction at large momentum transfer can be described perturbatively, whereas 
the quark-gluon structure of the initial and final hadrons still remains nonperturbative. Sometimes this physical idea 

^ ' can be put on solid theoretical ground. The so-called factorization theorems represent observable physical cross sections 
or amplitudes in terms of convolutions of perturbatively computable coefficient functions and nonperturbative objects 
like parton distributions, fragmentation functions, distribution amplitudes, generalized parton distributions, etc. 
The rich realm of factorization theorems can be divided in two different groups: 

■ 1) factorization of cross sections [usually associated with inclusive (semi- inclusive) processes], 
' 2) factorization of amplitudes in exclusive reactions. 

The border between the two groups is not absolute. The basic principles of quantum theory include the concept 
of the density matrix which allows one to combine the description of probabilities with the physics of interference. 
Another connection between the two classes of factorization theorems is provided by the optical theorem which gives 
us a relation between amplitudes and cross sections. 

A classical example of factorization for cross sections is deeply inelastic scattering where the required information 
on nonperturbative physics can be condensed into parton distributions. In the second group of hard processes with the 
factorization theorems for amplitudes an important role is played by the deeply virtual Compton scattering (DVCS) 
which can be described in terms of the generalized parton distributions (GPDs) 0,Qi0j0,|EIEI3>I1)0) OH El Ei3 • 
I ' In spite of the different physical meaning of usual parton distributions (probabilities) and GPDs (amplitudes), their 
Oh, theoretical description is based on the matrix elements of the same bilocal light-ray operators. For example, in the 
2^ ' case of quark distributions one deals with the operators 



o 

o 
in 



- 1—1 

X 



O r (A,n)=^(-^Wf^) > ( L1 ) 



2 J V 2 

where T is some spin matrix and the quark held ip is taken along the light-like vector n 



n 2 



= 0. (1.2) 



The GPDs are defined in terms of nondiagonal matrix elements of the bilocal operator 

J^ x (P 2 \0 T (\,n)\P 1 ) (GPDs), (1.3) 

whereas the definition of usual parton distributions is based on diagonal matrix elements of the same operator. These 
diagonal matrix elements can be obtained by taking the limit P 2 — > Pi in Eq. I|1.3fl . Therefore the usual parton 
distributions are sometimes called forward parton distributions (FPDs). Thus, FPDs can be expressed in terms of the 
matrix elements 

^e iAa: <P 1 |Or(A,n)|F 1 ) (FPDs). (1.4) 

In this paper we will concentrate on the so-called helicity- flip GPDs. These GPDs correspond to operators Or with 
matrices V changing the helicity of the parton. Helicity-flip GPDs were introduced by Hoodbhoy and Ji in Ref. [l4| . 
Their complete classification (for the twist-two case) w as g iven by Diehl in Ref. [l5( . Phenomenological applications 



of the helicity-flip GPDs are discussed in Refs. [iH 111 IToT ITl fl7L . 
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B. Positivity bounds on GPDs 



Now we want to describe the main idea standing behind the derivation of the positivity bounds for GPDs. 
precise classification of GPDs (which is discussed in Sec. II Gfl is not important for this brief derivation. 
Let us consider the states 



The 



\ xaPS) = J f/ xX ^ (t) lPS) ' 



where a and S stand for the quark and hadron polarizations. 
Using the positivity of the norm in the Hilbert space of states 



for the linear combinations 



!<=> = £< 



> o 



i &m Pm Srn ) 5 



(1.5) 



(1.6) 



(1.7) 



we find 
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This inequality can be rewritten in the form 



c h c 



(x k a k P k S k \x 



Pm. Sii 



< 



k^m 



k 



\ c k\ 



(x k a k P k S k \x k a k P k S k ) 



(1.8) 



(1.9) 



The nondiagonal matrix elements on the left-hand side (LHS) are associated with GPDs whereas the diagonal matrix 
elements of the right-hand side (RHS) correspond to FPDs. Using the freedom of choice of arbitrary coefficients c k , 
one can derive positivity bounds on GPDs from inequality ljl.9|l : 



|GPD| < / (FPD) 



(1.10) 



where function / depends on the set of coefficients Cj used in the derivation of the inequality. 

Various inequalities for GPDs have been derived in Refs. QSIIlEllIllillllllillH 
shown in Ref. j34j |. these inequalities can be considered as particular cases of a general positivity bound which has a 
relatively simple form in the impact parameter representation for GPDs [2H I2H IsoL Isil I32I Isflj . In the brief derivation 
of positivity bounds sketched above, we ignored the problems of the light-cone singularities, renormalization, gauge 
invariance, etc. These issues are discussed in Refs . pil l36 j . The positivity bounds are stable under the one-loop 
evolution of GPDs to higher normalization points |3 4 | . The positivity bound of Ref. [34| was explicitly checked for 
one-loop GPDs in various perturbative models [37]]. The solutions of the combined positivity and polynomiality 
constraints are studied in Refs. [37l l38|. The positivity bounds can be used for self-consistency checks of models of 

gpds m m EJ IH El El • 



C. Standard classification of twist-two quark GPDs 

The standard classification of the leading-twist quark FPDs includes three distributions 

FPDs: q,A L q,A T q (1.11) 

known as unpolarized, polarized and transversity, respectively. Any of these FPDs can be represented as the forward 
limit of some GPD but the inverse statement is not true. For some GPDs the forward limit cannot be associated 
with FPDs. Therefore the number of leading-twist GPDs is larger than the number of FPDs. The full set of the 
leading-twist quark GPDs includes eight functions |l5j 

GPDs: H q ,E q ,H q ,E q ,H^E^,H^,E^. (1.12) 

The twist-two operators correspond to the following constraint on the matrix T appearing in Eq. i|l.lf> 



r(n7) = («7)r = 



(1.13) 
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The solutions of this constraint are 

r = (717), (717)75, [(a 7 ), (r7 7 )] , (1.14) 

where vector a should obey the condition 

(an)=0. (1.15) 

The three cases ((1.1411 correspond to three usual FPDs (|1.11|) . 

The set of GPDs associated with the same light-ray operators is larger. The covariant decomposition of the nucleoli 
matrix elements of operators p. 1(1 leads to eight twist-two parton GPDs 1)1.12(1 . These GPDs are listed in Table [I] 



TABLE I: Twist-two quark GPDs and FPDs 



r 


GPDs 


FPDs 


(717) 


H q .E q 


<7 


(717)75 


H q ,E q 


A L q 


[(07), (77,7)] 







We assume the normalization of the light-like vector n 

n{P 1 +P 2 )=2 

and the standard notation 

A = P 2 - Pi , 



P=\{Pl+P2) 



Then the precise definitions of quark GPDs are 



ge^([/(P 2 )|^(-^) (777)^ (^lin/M; 



--U{P 2 ) H"(nj) + —EHa^n fl A u U{P X ) , 
2m 

^(U(P 2 M (~) («7)75^ (y) MPl)) 



U(P 2 ) 



2m 



J d ^{U{P 2 )\\i> [(a 7 ),(n 7 )]^ ( MP!)) 



U(P 2 ){-H^[(a 1 ),(nj)] + H' T 



2 J 

(Pn)(Aa) - (An) (Pa) 



7/7. 



^ q ( 7 n)(Aa) - (An)( 7 a) ~ q (jn)(Pa) - (Pn)( ia ) 1 

t ^ T ~ r^lPu- 



(1.16) 
(1.17) 

(1.18) 



(1.19) 



(1.20) 



(1.21) 



m J 

Here m is the nucleon mass, and a is an arbitrary vector obeying condition ((1.15|) . 

The GPDs depend on the standard Ji variables [24| x [which appears in the exponent e lXx on the LHS of Eqs 
and e, t: 

n(Pi - P 2 ) 



* n(Pi + P 2 ) 
At fixed £ the allowed region for the variable t is 



t = (P2 - PiY = A' 



where 



to 



t < t < , 
4to 2 £ 2 



(1.22) 
(1.23) 
(1.24) 
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D. Compact notation for kinematic variables 

Most of the positivity bounds for GPDs look rather cumbersome if one writes them directly in terms of the standard 
Ji variables x, £, t. Both the derivation and the final representation of the positivity bounds become much simpler if 
one introduces several new variables. 

The first pair of auxiliary variables 

is well known. The variable x\ (a^) has the meaning of the momentum fraction of the initial (final) quark with respect 
to the initial (final) nucleon. The positivity bounds are derived for the region 



where the variables Xk obey the constraint 



We also define the variables 



\£\<x<l, (1.26) 
0<x k <l. (1.27) 

(1.28) 



2m 1 -t 



(1.29) 



y/(t -t)(i- e) ato-t 

In terms of these variables we have 

t Q = -4a 2 m 2 £ 2 , (1.30) 

t = -4a 2 m 2 (iy- 2 + £ 2 ) , (1.31) 



to 



to-t 



_ ,2 t 2 



t -t 



= , (1-32) 



m 2 v 2 



t - t Aa 2 



E. Short guide through positivity bounds 



(1.33) 



The literature on the positivity bounds for GPDs is rather extensive. Without trying to give a complete review 
of the accumulated results, let us make several comments which may be helpful for the orientation in the realm of 
inequalities for GPDs. 

Starting from the general inequality ljl.8JI . one can reduce this inequality to others using various strategies: 

1) One can concentrate on the case when the linear combination (|1.7fl contains only two configurations for the 
quark-hadron kinematics (x±, P\ and X2, P2) but includes all possible quark and nucleon polarizations. This way is 
chosen in this paper. 

2) One can include all possible values of x m , P m in the linear combination l|1.7[) but concentrate on the minimal 
set of quark and nucleon polarizations. This leads to the positivity bounds in the impact parameter representation 
which were derived in Ref. |34[ for helicity-nonflip GPDs. 

3) One can choose an intermediate way by taking the parton-hadron configurations x, P with fixed x (but arbitrary 
P) and work with arbitrary polarizations. This leads to the positivity bounds for GPDs with £ = 0. The case of 
helicity-flip GPDs was considered in this context in Ref. |33l |. 
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II. EXAMPLES OF THE RESULTS 
A. Simple inequality for GPD H%. 

As was already mentioned, our aim is to derive bounds on helicity-flip GPDs which have the structure (|1.10(1 . In 
order to give a preliminary impression about the results obtained in this paper let us consider one example: 



\H q T (x,i,t)\ < Jq 



x + £\ _ fx — £ 



(i-m 



4m 2 ( 2 



-l 



(2.1) 



On the LHS we deal with the helicity-flip GPD (|1.21[) depending on the Ji variables x, £,t ^1.22(1 . On the RHS 
we have the product of two unpolarized quark FPDs q taken at values x k (|1.25|) : 

qk = q{x k ) . (2.2) 

This positivity bound (as well as other positivity bounds) can be used in the domain i|1.26[l . where the variables Xk 
belong to the physical region < x k < 1 l|1.27|l and the FPDs q(xk) do not vanish. In terms of variables a (|1.28|) . v 
(Q§|) and q k (j22l we can rewrite inequality l|2.1|l in the compact form 



\H q T \<a^{l + ^e) qi q 2 . (2.3) 
Note that this positivity bound for the GPD Hj, contains only the unpolarized FPD q. 



B. Enhanced positivity bound on H\ 

In fact, using the polarized FPD A^q in addition to the unpolarized FPD q, one can derive another bound for the 
same GPD: 

\H q T \ < f [(rj> 2 - + r+4) 2 + v^ 2 {r+r- + r^f] ^ . (2.4) 
Here we use the compact notation 



r± = yjq{x k ) ± A L q{x k ) . (2.5) 

Note that the square root on the RHS is applied to a positive expression since 

\A L q\<q. (2.6) 

Let us show that the previous inequality l|2.3|l is a trivial consequence of the stronger inequality (|2.4(1 . Indeed, 
starting from the general inequality 



we find 



+ a 2 b 2 \ = |a- b| < |a||b| = ( ai ) 2 + (a 2 ) 2 y/ (h) 2 + (b 2 ) 2 , (2.7) 

r^r-+r+r+ < ^ {r^f + {r+f ^ {r^f + {r+)\ (2.8) 

r x +r 2 - + r x "r+ < (r^) 2 + (r+) 2 ^(r") 2 + (r+) 2 . (2.9) 

According to Eq. I|2.5|l we have 



r,. 2 + (r+) 2 = 2q k . (2.10) 



' k ) ^ V k 

Taking a linear combination of inequalities Q2.8[l . (|2.9|l and using Eq. I|2.10|l . we find 



(rrr 2 -+r+r+) + v 2 ? (r+r^ + r^r+) < 4 (l + ^ 2 £ 2 ) q x q 2 . (2.11) 
Using this inequality, we see that positivity bound (|2.3|l is a consequence of inequality (|2.4I) . 
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C. Positivity bound on using transversity FPD 



Inequality l|2.4[) can be also enhanced. In order to write the enhanced positivity bound in a simple form we introduce 
a special notation. Starting from FPDs Ijl.lip . let us construct the following combinations 



(2.12) 



( Q x \ 




/V 


-A L q + 2A T q \ 


Q 2 






q- a l <7 


Q 3 




q- 


- A L q - 2A T q 


\Q 4 J 




{ 


q-^ L q J 



One entry is duplicated here: 



but this redundancy will be helpful for our later work in Sec. Mil 
Using Eq. I|2.6|l and the Soffer inequality jl^ 

2|A T g| <q + A L q, 

M 



we conclude that all functions Q are positive 



It is convenient to introduce the notation 



with Xk defined by Eq. Ijl.25jl . 



Q M > . 



nM 
Hi. 



^Q M {x k ) 



(2.13) 

(2.14) 
(2.15) 
(2.16) 



Then our inequality for H!L has the form 



- 4 

fc=l,3 



[(J2*JJ* 



R1R2) 



u 2 e 



(-n, 1 K 2 + K 1 K 2 ) 



1/2 



(2.17) 



In Appendix we show that the previous positivity bound l|2.4|) is an algebraic consequence of inequality (|2.17|) . 

Let us summarize. We have presented three inequalities (|2.1|) . I|2.4|l and 12.17|l for H T . Inequality l|2.17|l is the 
strongest one, it can be used as a starting point for the derivation of inequality (|2.4|l . In turn, positivity bound 12.11) 
follows from Eq. (|2.4|) . Thus, the problem reduces to the derivation of the strongest bound l|2.17|l . We postpone this 
derivation till Sec. ITvTI 

Although formally the bound l|2.17|l is the strongest one, in practical applications it is not too helpful since it 
contains the transversity distribution A^q which is poorly known. In Appendix iBl we show how the maximization of 
the RHS of the bound (|2.17l) with respect to the unknown transversity distribution Apq leads to the weaker inequality 
(|2.4() . This weaker bound contains the forward distributions q and A^q on the RHS. Depending on the availability of 
phcnomenological data on Al<? one can either use this bound directly or one can maximize its RHS with respect to 
A^q. In the second case one arrives at the weakest bound (|2.1|) for H^. 

We have described the hierarchy of positivity bounds the GPD Hj,. The case of other GPDs is similar. Below we 
will derive the bounds of the strongest type (|2.17|l for all helicity-flip GPDs. After that we will show how weaker 
bounds can be derived by eliminating Apq and A^q. 



III. GENERAL POSITIVITY BOUNDS FOR GPDS 



A. Quark case 

In Sec. |n]we have considered several examples of positivity bounds. Now we want to turn to the derivation of the 
inequalities. The idea of this derivation is represented by Eqs. I|1.8ll i|1.10fl . As for the technical realization of this 
idea, an essential part of the work has already been done in Ref. [28| , where a general positivity bound was derived 
for the so-called helicity amplitudes. In this section we briefly describe the main result of Ref. psf. In Sec. IIVI we 
will turn to the analysis of the special case of positivity bounds for the helicity-flip GPDs. 
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The general positivity bounds on GPDs were derived in Ref. [28j using a special matrix notation for GPDs and 
FPDs: 

1) matrices F q (x) made of FPDs q, A^q, Arq, 

2) matrices Bj(x,£,t) made of GPDs H*,E«,H*,E*,H%,E«.,H*,E* JET2J. 

Both F q and B q are 2x2 matrices. In addition to the matrix structure, F q and B% have the subscript s which 
runs over the values s = 1, 2. The precise definition of F q and B q s will be discussed later. 

Using these GPD matrices B q (x,(;,t) and FPD matrices F q (x), one can write the general positivity bound derived 
in Ref. [H in the form 



<^Tr([F«( Xl )L s F«(x 2 )Lt] 1/2 ) . (3.1) 



This inequality holds for arbitrary 2x2 matrices L s (s = 1,2). L\ stands for the matrix Hermitian conjugate to L s . 
The power 1/2 on the RHS of inequality l|3.1|) should be understood in the matrix sense. The variables Xk (k = 1, 2) 
are given by Eq. i|1.25fl • 

Using the freedom of choice of matrices L s , one can derive bounds for arbitrary GPDs from the inequality (3.1). In 
Ref. [2^ the general positivity bound (|3.1(l was used for the derivation of inequalities for the GPDs H q , E q , H q , E q . 
In this paper we will derive inequalities for the helicity-flip GPDs Hj, , Ej, , Hj, , Ej, from the same general inequality 

Now we turn to the definition of matrices F q and B q which appear in the general positivity bound (|3.1|l . Matrices 
F q (x) are defined via the linear combinations Q M (x) of FPDs l|2.12|l : 



( Q1{X) 9 ° I , F«(x)=( 4 ° ] . (3.2) 

1 Q 2 {x) J ' 2V ' \ Q\x) 1 



The definition of matrices B q (x, £,t) can be found in Appendix IA 61 In this section we do not need the explicit 
expression for B q (x : ^ : t). But it is important to understand that the two matrices B q (with s — 1,2) contain 
altogether 4 + 4 = 8 linearly independent matrix elements. This is exactly the number of the full set of twist-two 
GPDs (|1.12(l . Hence the matrix elements B l J of matrices B s can be used as an alternative basis in the space of linear 
combinations of GPDs. Thus, any GPD (and any linear combination of GPDs) G can be represented in the form 

2 

G q (x,U) = E L{\G q )(B q s r (x,U)- (3.3) 

i,j,s=l 

where L{ l (G q ) are coefficients of the decomposition. Interpreting the coefficients L l J(G q ) as matrix elements of 2 x 2 
matrices L s (G q ) we can write 

2 

G*(x,t,t) = X>[£.(GW(*>&*)] • (3-4) 

a-l 

Inserting this decomposition into the LHS of inequality l|3.1|l . we find 

2 

\G q (x,U)\ <^{[^(^) L s(G v )FK^)Li(G q )} 1/2 } . (3.5) 

3 = 1 

In order to compute the RHS of inequality (|3.1|l . we use the general formula 



Tr (5 1/2 ) = V TrS + 2Vdet S (3.6) 
valid for 2x2 matrices S. It follows from the fact that any matrix T obeys its characteristic equation 

P(A) =det(T- A), P(T) = 0. (3.7) 
In the case of 2 x 2 matrices, this results in 

det T - (TrT)T + T 2 = . (3.8) 
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Taking the trace of this equation and substituting T — S* 1 / 2 , we obtain identity l|3.6|) . Using identity l|3.6[) and short 
notation R^ 1 H2.16[) . we easily compute the RHS of inequality (|3.5|) : 



1/2 



\G«\<\ ]T (Rr i R2) 2 \[Li(G^} M J + 2\detL 1 (G^)\RlR 1 2 RlRt 



M,N=1 
2 



1/2 



Y, « +2 < +2 ) p 2 (G<0] MJV | 2 + 2 |detL 2 (G*)| i^i?^ 



M,N=1 



B. Gluon case 



(3.9) 



The generalization for the case of gluon GPDs is straightforward. We follow the definition of gluon GPDs used in 
Ref. 15]. In the gluon case, inequality 1)3.1(1 should be modified by a factor of \J\ — £ 2 (see Ref. [H)): 



z 



/2 



s=l 



The gluon analogs F^(x) of quark matrices (|3.2[1 are 



where we use notation 



Ff(x) =Fi(x) = I X9+ ^ ° 

V xg (x) 



g ± {x)=g(x)±Ag(x). 



(3.10) 



(3.11) 



(3.12) 



Here g(x) and Ag(x) stand for the forward gluon unpolarized and polarized distributions, respectively. Note that in 
the case of spin-1/2 hadrons we have no forward gluon transversity distribution. Thus, the gluon analog of inequality 
(|3.9|) can be obtained by replacing 



\G«\ 



1 



\G 9 \ 



on the LHS of inequality (|3.9|) and replacing 



(Rl) 




(vt\ 


Rl 
Rl 




Vk 

vt 


\RtJ 




\Vk 1 



(3.13) 



(3.14) 



on the RHS. As a result, we obtain 



V 



1 5 s=l I M,JV=± 



1/2 



J MAT I 



2|detL s (GS)|y+y 2 V% 



Here we use the values M, N = ± for matrix indices instead of 1,2. 
According to Eq. i|1.25f) we have 



Therefore 



v^vi^e = Vx T ^e 



\GV\<V^eY\ E (y¥y%) 2 \[Ls{GV)] MN \ 2 + 2\AetL s {&)\y+yly^ y2 - 
s=l I M,iV=± 



1/2 



(3.15) 



(3.16) 



(3.17) 



(3.18) 
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IV. DERIVATION OF POSITIVITY BOUNDS FOR HELICITY-FLIP GPDS 

In order to apply the general positivity bounds l|3.9|) and l|3.18[l to some GPD G, we must know the corresponding 
matrices L S {G). In Appendix I A 71 we describe how L S (G) can be computed using Eq. I|3.4[l . Below we list the results. 



A. Matrices L S (G) for quarks 



Using notation 



we can write 



Vs = 



4, s = 1, 
-1, s = 2, 



4 ^ -u£ -ri a J 



(4.1) 



(4.2) 



L s [H* 



t£_ I 

4a 1 f] s 



(4.3) 



L s (E« ) = — 



-l + £ 



4 V 1 + C -2%^ / ' 



(4.4) 



av I l-£ 
4 I 1 + £ -2r?X 



(4.5) 



B. Matrices L S (G) for gluons 

In Appendix I A 71 we show that matrices L S (G) for gluon GPDs H^, H^, Ej,, Ej, differ from the corresponding quark 
matrices by the factor of v. 



Therefore 



L S (G 9 ) = vL s (G <1 ) (G 9 = H 9 T ,H 9 T ,E 9 Tl E 9 T ). 

av ( r) a -v£ \ 



L s {H 9 T ) = 



-H -r} s 



- i ( 



V s I 

Vs 



(4.6) 



(4.7) 



(4.8) 



L. (E 9 T ) 



9\ - 



w 2 I + £ 

4 I 1 + e -2ry s2 v 



(4.9) 



I o 1 -£ 

4 I 1 + £ -2t?X 



(4.10) 



Note that in the positivity bound (|3.18|) we use the values ± for the matrix indices of L s . The above explicit matrix 
expressions assume the interpretation "+ — 1" , "— = 2". 
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C. Example: derivation of inequalities for H^, 



Let us show how the general inequality (|3.9|l can be used for the derivation of the positivity bound for the GPD 

I) into Eq. we find 



Hip. Inserting the explicit expression (|4.2|l for L S (H, 



W q T \ < \ 



1/2 



(rIr 2 



According to Eqs. I|2.13[l and (|2.16() we have 



R 1 R 2 ) 



1/2 



(4.11) 



i?f — i? 2 . 



(4.12) 



Now we see that inequality l|4.11|l coincides with the inequality (|2.17|l which was already discussed above. As was 
explained in Sec. [HI two other positivity bounds for Hi, l|2.3[l and l|2.4|l can be derived from the already proved 
inequality l|2.17|l . 



V. RESULTS 



A. Hierarchy of positivity bounds 

In this paper we have chosen quark GPD Hj. in order to illustrate the methods used for the derivation of positivity 
bounds for GPDs. We have derived three inequalities: 

1) "strong" inequality (|2.17l) where the GPD is bounded by a combination of all FPDs (transversity, longitudinally 
polarized and unpolarized) ; 

2) "intermediate" positivity bound l|2.4|l where the GPD is bounded by a combination of unpolarized and longitu- 
dinally polarized FPDs but the transversity FPD does not appear; 

3) "weak" positivity bound l|2.3|l which involves only the unpolarized FPD. 

From the mathematical point of view we have a set of inequalities which can be derived one from another, starting 
from the strong inequality l|2.17f) . However, currently the use of the strong positivity bound is limited by the absence 
of reliable data on the transversity FPD. On the other hand, the intermediate and weak positivity bounds do not 
contain the transversity FPD and can be used in practical applications. 

The above classification of positivity bounds for GPDs was described in terms of the quark GPD H!^ but the 
generalization for other quark GPDs is straightforward. 

In the case of gluon GPDs the situation is simpler. Since there is no gluon transversity FPD for the nucleon, we 
have only two types of positivity bounds for gluon GPDs: 

i) "strong" positivity bounds containing both polarized and unpolarized gluon FPDs, 

ii) "weak" positivity bounds involving only the unpolarized gluon FPD. 
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B. Strong positivity bounds for GPDs 

1. Strong positivity bounds for quark GPDs 

As was explained above, we use the word "strong" for those inequalities which come directly from the general 
positivity bound (|3.1|) and contain all types of FPDs on the RHS. The set of "strong" inequalities for quark GPDs is 



\H%\ < 



k R\ + R\Rl) 2 + v 2 e (R\Rl + RlR k 2 ) 2 



fe=l,3 



1/2 



\Hl\<—R\Rl, 



fc=l,3 



fe=l,3 L 



d/c d2 d2 r?/c 



4^ 2 a 4 (i? 2 i? 2 2 ) 



1/2 



4^ 2 a 4 {RlRlf 



1/2 



(5.1) 
(5.2) 
(5.3) 

(5.4) 



Parameters a and i/ are given by Eqs. Qgj) and fTSity . The first inequality iJOJ is nothing else but Eq. lj2~T7|) 
derived above. The derivation of the other inequalities is similar. In all these inequalities, the GPDs on the LHS are 
taken at values x, £, t [e.g. Hj, = Hj,(x,£,i)}. The FPDs appear in these bounds at points x k (|1.25|) via the quantities 
R\. According to Eqs. I|2.16[) and (|2.12|) we have 







Rl 


-( 


\Rl) 





(5.5) 



y/g{xk) + A L q(x k ) + 2A T g(x k ) 

\Jq{xk) - A L q(xk) 
yjq{x k ) + A L q(x k ) - 2A T q(x k ) 

2. Strong positivity bounds for gluon GPDs 

The difference between the quark and gluon cases is so small that the inequalities for gluon GPDs can be directly 
obtained from quark inequalities l|5.1l) - (|5.4l) by a simple replacement. The modification consists of two factors: 
1) factor 



(5.6) 



arising from the difference between the underlying general quark l|3.9|) and gluon i|3.18[l inequalities, 
2) factor v coming from the relation (|4.6(l between the quark and gluon matrices Lk(G). 

Taken together, these two factors give vw. Now we must insert this factor into the quark inequalities 15.1(1 (|5.4|) 
and replace 



(5.7) 



(Rl\ 


(vt\ 


Rl - 


y k 


\Rl) 





where according to Eqs. (|3.12() and (|3.15|) we have 



Vk = Vdixk) ± Ag(xk) ■ 



(5.8) 
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As a result, we obtain 



\m < ^ 



(vt vt + Vi V2 Y + ^ 2 £ 2 {vi vt + viv2 Y 



1/2 



\m\ < 



2a 



1^1 < ^ 



2a 



1 Tl - 2a 



Vi 2/2 , Vi Vi 



+ 



yi 2/2 , Vx 2/ 2 



1 1/2 



+ 4^a 4 (yry 2 -) 5 



+ 



+ 4, 2 ar(2/r%) 2 



1/2 



(5.9) 
(5.10) 

(5.11) 
(5.12) 



C. Consequences from strong positivity bounds 



1. Positivity bounds for quark GPDs without forward transversity distribution 



The elimination of the transversity FPD from the positivity bounds for quark GPDs is discussed in Appendix IB1 
Here we just list the results: 



H q T \ 


< 


-[( 

2 V 


r i r 2 + r i 


r 2 


) 2 






v 2 
— r 
2a 








H q T \ 


< 


1 r 2 ) 






\E q T \ 


< 


V 


( r i r ? 


r l 


'2 


2a 


[{iA + 


1 - 


-s 


\E q T \ 


< 


V 


( T l r 2 


r l 


'2 


2a~ 




1 - 


-i 



-\2 



1/2 



+ Av 2 a A (r x r 2 )' 



1/2 



+ A V 2 a^ 2 (r-r-y 



1/2 



(5.13) 
(5.14) 

(5.15) 
(5.16) 



Parameters rjf appearing on the RHS were defined in Eq. i|2.5|l . They depend on FPDs q and A^q only. 



2. Positivity bounds for GPDs containing only unpolarized FPDs 



The positivity bound for Hj, l|2.3(l containing only the unpolarized FPD has been already derived from Eq. (|5.13(l 
in Section lilBI A similar inequality for can be obtained from inequality Eq. (|5.14l) by noticing that < \J2q k . 
Thus, for the quark GPDs we have 



with the qt defined in Eq. I|2.2|l . 
For gluons we similarly find 



where g k = g{x k ). 



\H*.\<ay/(l + v 2 ?) qi q 2 , (5.17) 
\H«\ < v 2 a~ 1 ^T 2 (5.18) 



\H 9 T \ < vwa^{l + v 2 e)9i92, (5.19) 
\H 9 T \ < v^a^Wy/gTgl, (5.20) 



VI. OTHER BOUNDS AND APPLICATION TO DVCS 



Helicity-flip gluon GPDs app ear in the QCD description of DVCS (via the interference of the DVCS amplitude with 
the Bethe-Heitler process) [l(J,[l5|. Helicity-flip quark GPDs can also be accessed in the diffractive electroproduction 
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of two vector mesons 0,^3^3- I n practical applications one usually deals with linear combinations of GPDs. In 
principle, we can simply combine the trivial inequality 



C\Ht + CiHt + C^Et + C4-&T 



< \cxHj 



c 2 H 7 



\c 3 E T \ + 



c^Et 



(6.1) 



with the above bounds for separate GPDs H x , Et , Ht > Ex ■ However, one can derive a stronger bound for the same 
linear combination using the basic inequalities (|3.9[l . (|3.18|l . 

In order to illustrate this method, let us consider the interference of the Bethc-Hcitler and DVCS amplitudes in 
the case of unpolarized protons and electrons. The cos3</> interference term can be expressed via the following linear 
combination of gluon helicity-fhp GPDs [TnL ITHI | 



,-DVCS 



F 2 H T - F X E T - 2 LFi 



Am? 



f 2 m 



(6.2) 



Here Fi t2 are Dirac form factors. In order to apply the general positivity bound (|3.18|l to G^ cs , we first must 
compute the corresponding matrix L S (G^ )'■ 



L S {G 



DVCS^ 
3<p ) 



F 2 L S (H T ) - FiL s (E 9 T ) -2 [Fx 



Am 2 



F 2 L„ H, 



Inserting Eqs. (|4.7I) - (|4.9() and using expression l|1.31|l for t, we find 



(6.3) 



L s \S"iA> ) — —r 



v K*2 + (1 - F x \ 
4 I -v \iF 2 + (1 + Fx] Vs [(l + 2v 2 e)F 2 + 2F 1 v 2 e] 



Let us apply the general formula l|3.18|l : 



1/2 



M„.N\ 



s=l I M,N=± 



MN 



2|detL s (G^ CS ) 



vi vi Vi y 2 



Using notation 



gf = g*(x k ) = (jfi) 2 , 



we obtain 



E ( 

M,N=± 



MN 



^) 2 { {F2?gigt + [-v£F 2 + (1 - OFiv] 2 gtg 2 + [-v£F 3 - (1 + H) Fiv} 2 gt g^ 
[2F lV 2 e+F 2 {l + 2v 2 e)] 2 9i9 2 } , 



(6.4) 



(6.5) 



(6.6) 



(6.7) 



det L s (G 3rj) 
As a result, we find from Eq. I|6.5|l 



(F lV a^) 2 + {F 2 f(l + v 2 e) 



(6.8) 



| r DVCS| < WaV { ( F \ 2 n + n + 

\G 3lp I < — — I \*2) 5i 9 2 

+ v 2 m - (1 - F,f g+g- + v 2 [£F a + (1 + F{\ 2 g+g- 

+ [2F lV 2 e + F 2 (i + 2 v 2 e)] 2 g 1 g 2 



{F lV a- l ) 2 + {F 2 f(l + v 2 e) 



1/2 



(6.9) 
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VII. CONCLUSIONS 



In this paper we have derived a set of positivity bounds for the quark and gluon helicity-flip GPDs. For every GPD 
we have several positivity bounds ranging from the simplest bounds in terms of unpolarized forward distributions to 
the stronger bounds involving the transversity forward distribution. 

It should be stressed that the method developed in this paper can be more useful than the bounds presented here. 
Indeed, in practical applications one usually deals with linear combinations of GPDs depending on the specific process. 
The method described in this paper allows one to derive positivity bounds for any linear combination of GPDs. 

The power of positivity bounds is limited to the region |x| > £. Although we have no explicit positivity bounds at 
|x| < £, this region is not completely isolated and free of any constraints. Indeed, the polynomiality of GPDs imposes 
restrictions involving both regions. It should be also emphasized that the inequalities derived in this paper are just 
a small part of the general positivity bounds which can be written in terms of the impact parameter representation 

M 

Thus, one has an infinite amount of the polynomiality and positivity constraints (in the impact parameter repre- 
sentation). The construction of the complete general representation for GPDs obeying all these constraints remains 
a challenging problem. 
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APPENDIX A: HELICITY AMPLITUDES 



1. Definition 

It is convenient to analyze the general positivity bounds for GPDs using "helicity amplitudes" introduced in Ref. [l5| 



A' fi' ,A/x 



dz 



(p',X'\0^^(z)\p, A) 



(Al) 



!+=0„ 



Here M '. M are bilinear quark light-ray operators with the polarization indices u, fjf. Next, \p, A) is a nucleon state 
with momentum p and polarization A (in the sense of light-cone helicity states |46(), P + is the light-cone component 
of vector P = (p + p')/2. 

With two values for each polarization index of A\i „/ xu we have 2 4 = 16 components but due to the parity invariance 

M 



A_A',-/i';-A,-jU — ( — 1) ^ "^Ax'/i'jAju 

only eight components are independent. This is exactly the number of twist-two GPDs l|1.12fl . 



(A2) 



2. Quark helicity amplitudes 



We take the expressions for the helicity amplitudes A\, , Xfl via GPDs from Ref. 0] using variables a (|1.28|) and 

1 



Al 



A%.- 



A'L 



2a 
1 

2a 
a 

~2^ 



H q + H q \- a V [E q + E q 



E q - E q 



A q 



H q - H q 
(E q - (,E q 



* 2 e ( 



(A3) 
(A4) 
(A5) 
(A6) 
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At. 



A q 



Hi 



Hi 



2 



(e*. + E q T ) 



E q T -E q T 



+■+- 



a (a- 2 H<<, + V - 2 Hl - i 2 E q T + £Ef, ) . 
-av~ 2 Hl . 



(A7) 

(A8) 

(A9) 
(A10) 



3. Gluon helicity amplitudes 

Now we list the expressions for the gluon helicity amplitudes Ay , ^ via GPDs from Ref. |Tf| (taking e = +1 in 
equations of Ref. 15]). Amplitudes which conserve the gluon helicity have the same form as in the quark case: 

(All) 
(A12) 
(A13) 
(A14) 



^++,++ 


1 - 

~~ 2^ . 


^-+.- + 


1 " 


^++.-+ 


a 

~ ~2^ 




= -( 



H 9 + H 9 j - a 2 i 2 [E 9 + E 9 
H 9 - H 9 ) - a 2 i 2 (e 9 - E 9 
' E 9 - $E l 



but the helicity-flip amplitudes differ by a factor of v 1 

A ++,+ 
A'i 



H 9 T 



V (El + K 



-+,- 



H 9 T + i-^ ( E", E" 



A 9 ++ ^__=av- 1 [a- 2 H 9 T 



v- 2 H 9 T ~t 2 E 9 T +tE 9 T] . 



av- z H 9 T . 



(A15) 

(A16) 

(A17) 
(A18) 



4. Expressions for quark GPDs via helicity amplitudes 

Solving the system of Eqs. (|A7|) - (|A10p with respect to GPDs, we find 



Hi 
Hi 
El 



(A q 



- x v 2 A q . 



2av 2 A q 



-+,+- ' 

av[(l-OAl +i __ + (l + 0^+,+-J ' 

E q T = av [(1 + - (1 " A -+,—] + ^v 2 A q _ +t+ _ 



(A19) 
(A20) 
(A21) 
(A22) 



5. Expressions for gluon GPDs via helicity amplitudes 

Now we want to solve equations <|A15(1 - <|A18|) with respect to the gluon GPDs. Note that the only difference form 
the quark equations is the extra factor of v~ l in Eqs. HA15() - (|A18() . Therefore we can read the gluon results form 
quark expressions (|A19|) — (|A22fl : 

H 9 T = av {A 9 ++ __ + A 9 _ + ^ + _) + av 2 £ (A 3 _ +> __ - A 9 ++>+ _) , (A23) 

Hi = -a- 1 v 3 A 9 _ + + _ , (A24) 

E 9 T = av 2 [(1 - Cs A -+,— + (1 + 0^++,+-] + 2av 3 A g _ +i+ _ , (A25) 

El = av 2 [(1 + €)A B ++r+ _ - (1 - A -+,—] + ^u 3 A 9 _ +>+ _ . (A26) 
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6. From A to B 



Instead of the tensor Ay/*',A/t) it is convenient to use the matrix A a b defined by the linear transformation 

A ab = ]T (V^')* A XI ^ M V^ (a,6=l,2,3,4), (A27) 

where has following nonzero components 

Vt ++ = V 3 ++ = V~ = -Vf~ = V 2 + ~ = = V± + = -V 2 + = 1 . (A28) 
Matrix A a \, has a block diagonal structure 



where the 2x2 matrices Si, B 2 are 



Now we find from Eqs. IA2I>. (|A27l 





A = 


( ill 


Il2 N 


~ \A 21 


^4.22 y 


.28)) and I 




A++,- 


i 

~ 4 


A-+,- 




A++,+- 


1 

~ 4 


A-+,+- 


1 

~ 4 



Si 

s 2 



(A29) 



[(Si)n - (S 2 ) u ] , (A31) 
| [(Si) 21 + (S 2 ) 2 i] , (A32) 
[(Si)i2 + (S 2 ) 12 ] , (A33) 



[(S 2 ) 22 - (S x ) 22 ] . (A34) 



7. Example of a calculation of Lk (G) 

Let us illustrate the calculation of matrices Lk (G) for the GPD Hj.. Inserting Eqs. i|A31fl - i|A34|) into expression 
(|A19|) for H%, we obtain 

H q T = \ {[(Si)u - (S 2 )n] + [(S 2 ) 22 - (Si) 22 ]} 
- Ja^{[(Si) 21 + (S 2 ) 2 i] + [(Si) i2 + (S 2 )i 2 ]} 



a 

4 Tr 



1 -^)b 1+ [ 1 ~* )b 2 



(A35) 



Comparing this result with the general representation (|3.4|l . we arrive at the result l|4.2|l for Lk {H^). The other quark 
matrices Lk (G q ) (|4.3fl - l|4.5[l can be computed in the same way using equations (|A20(1 ~ I|A22(1 for Hj,, E!^, Ej,. 

The gluon matrices Lk(G 9 ) can be immediately read from the quark matrices Lj~{G q ), using relation (|4.6|l . This 
relation follows from the fact that expressions l|A23|) - (|A26|) for gluon GPDs differ from the corresponding quark 
expressions (|A19|) - (|A22|) only by a factor of v . 

APPENDIX B: EXCLUSION OF THE FORWARD TRANSVERSITY DISTRIBUTION FROM THE 

POSITIVITY BOUNDS 

1. Bound on H\ 



In section III CI it was announced that inequality l|2.4|) can be derived from the stronger inequality (|2.17|) . In this 
appendix we will demonstrate this. 



Introducing the variables 



and the function 



h(v!,v 2 )= (l + V^f + v 2 t 2 (V^+V^) 2 



1/2 



we can rewrite inequality H2.17|l in the form 

\H*\ < jR\R 2 2 [h( Cll ,c 12 ) + h(c 31 ,c 32 )} . 
Obviously Cjk > here. Now we want to show that function h is convex: 



d 2 h(v 1 ,v 2 ) 



dv k dvi 



> 0. 



Let us define 



so that 



(vi,v 2 ) = (1 + \/viv 2 ) 2 + (v^) 2 (y^T+ Vv2) 2 



h(vi,v 2 ) = y/4>(vi,v 2 ) . 



Then 



d 2 h(vi,v 2 ) _ 1 



dvkdvi 4</> 3 / 2 



d(j)_d^_ _ d 2 q 



dvk dvi dv k dvi 



Thus, in order to prove the convexity of h, we must check that matrix 



is positive definite: 

The positivity of the first term is obvious: 
Now let us show that the second term 

is also positive definite. Since 
it is enough to show that 

Indeed, 



deb d 2 cp 

Oki = -^—t, 20- — — 

di'k ovi ovkOvi 



\D U \\ >0. 



dv k dvi 



> 0. 



-2(f> 



d 2 



dv k dvi 



> 



0>O, 



o 2 



dv k dvi 



> 0. 



0(wi,«a) = i + 2 l + K) 2 y/viVi + K) 2 (vi + v 2 
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so that 



dvkdvi 



1 r 
2 



y/VlV 2 



V2_ _^ 

V2 



(B15) 



Obviously this matrix is degenerate and non- negative. Combining this positivity with Eq. 1|B7|1 . we complete the 
derivation of Eq. (|B4|) . This result means that function h is convex. Therefore 



h{cn, ci 2 ) + /i(c 3 i, c 32 ) < 2h 
Together with inequality i|B3|l this yields 



H q T \ < -R\R\ [h( Cll ,c 12 ) + h(c 31 ,c 32 )} < -R\R\h 



Cll + C31 C12 + c 32 



01 D2 D 2. ( C H + C 31 C 12 +C 32 



According to Eqs. JbTJ), and lj2~T2l) we have 



Clfc + c 3fc _ 1 
2 ~~ 2 



*2 



*2 



(g + A £ g) A 



Therefore inequality (|B 1 T|> takes the form 



\H q T \ < ^r^h 



(B16) 



(B17) 



(B18) 



(B19) 



According to Eq. (|B2(1 this leads to inequality l|2.4|) . 



2. Intermediate bounds on and 

Now we want to repeat the work on the elimination of Aj^q for other quark positivity bounds. Note that the 
positivity bound H5.2J) for Hj, does not contain the transversity distribution. Indeed, the RHS of this inequality 
depends on FPDs only via i? 2 and according to Eq. (|5.5|l R\ is independent of At9- Thus, the problem of the 
elimination of the transversity must be solved only for GPDs Ej, and E\ . 

We can rewrite inequality l|5.3|) in the form 



\E q T \ < ^ 



v^7 



b 2 



(B20) 



where 



ak = — — — - H 

2a \ 1 + £ ^ 1 - £ 



(B21) 



b = v 2 aR\R\ . 



(B22) 



Using the general inequality 



(ai) 2 +& 2 + VM 2 +6 2 



< J (ai)2 + (fl3)2 +6 2 



(B23) 



we find from inequality (|B20I 



im< /( ai ) 2 + (a 3 ) 2 



6 2 . 



(B24) 
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Here 



(aif + (a 3 f = ^ 



R\R\ 



V 

4q,2 



^) 2 + (i??) 2 



(i-O s 



op 2 02 ^1-^2 + -Rj-Rj 
2^2^ j— ^ 



(B25) 



According to Eq. (|2.12(l we have 



1 r 



(Rl) 2 + (^) 2 = q(x k ) + A L q(x k ) = (r+) : 



(i? 2 ) 2 = q{xk ) - A L q(x k ) = (r-y 



(B26) 
(B27) 



R\R\ + R\Rl < J {R\Y + (RlY (R\y + (R 



2r+r+ 



Using these expressions, we find from Eqs. I|B25|I and l|B22 

(ai) 2 + (a 3 ) 2 <^ 



1 ' 2 



'2 ' 1 



(B28) 



(B29) 



b = v 2 ar^r^ . (B30) 



Now we insert these results into Eq. (|B24|) and obtain the positivity bound (|5.15|) . 

Note that the bound l|5.4|l on Elj, differs from the bound (|5.3() on Ej, by the replacement 



v 2 a 4 {RlRlf - e» 2 a 4 {R 2 R 2 ) 2 ■ (B31) 



Making this change in Eq. I|5.15|l . we arrive at the positivity bound l|5.16|) . 
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